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^^ . Abstract 



^SJ ' We calculate the damping rate 7; for longitudinal gluons with zero monien- 

Cn , turn in finite high temperature QCD and show that some of its contributing terms 

are infrared divergent. This is in contrast with the expectation that this damp- 
ing rate is to be equal to the corresponding one 74 for transverse gluons which 
is known to be finite. Our calculation was motivated by the fact that similar 



divergent terms occur when we calculated in a previous work 7^ to order p , p 
being the momentum of the gluon. After we present our results, we briefly discuss 
Oh, them. 



In a previous work [1], we undertook the calculation of the damping rate 7t(p) for transverse gluons 
with soft momentum p in finite- high-temperature QCD. Knowing that 74 at zero monicntum is finite [2], 
our main aim in that work was to argue that for small enough p = |p|, the expansion of 7t(p) in powers 
of [p/rag)^ is valid, rag being the inverse gluonic correlation length. We suggested that, when contrasted 
with a previously recent calculation [3] based on letting p come from the hard limit (set by the temperature 
T) down towards the interior of the soft region (set by mg ~ gT where g is the QCD coupling constant), a 
calculation that produced a behavior of the form ln(l/g) for jt{p), such a result could possibly be indicative 
of a kind of transition around some 'critical' scale. We therefore expanded jt{p) in the following manner: 



7t(W = ^T7^ '''0 + ati[—) +... , (1) 



247r 



and gave an analytic expression for an (the first coefficient ato was already calcualted in [2]). Nc is the 
number of colors. 

The natural step forward was to calculate numerically an- But quite interestingly, when we started 
manipulating the expression of ati we obtained, we realized that some of its contributing terms are infrared 
divergent (in a sense that will become clear later in the text). Thus, on the one hand, the quick conclusion 
is that ati may be infinite and hence, the expansion in (1) is not valid after all. On the other hand, such 
a potential infrared divergence seems in contradiction with the physical expectation that the quark-gluon 
plasma is stable for, at least, very small gluonic momenta. 

In order to look into the matter more closely, we undertake in this work the derivation of 7; (0) , the 
damping rate for longitudinal gluons at zero momentum. Our motivation is that, first, this damping rate 
is to be equal to 7t(0), simply because at zero momentum, there is no difference between transverse and 
longitudinal directions. Second, the interesting point in this calculation is that in order to get 7; at zero p, 
it is already necessary to expand it to order {p/irig)'^. After the calculation is performed, we indeed recover 
the finite analytic expression that is equal to 7t(0), but, surprisingly enough, together with additional terms 
quite similar in form to the ones we obtained in the case of 74 (p) , terms when taken individually present an 
infrared divergence. To the best of our knowledge, in this context, such infrared-divergent terms have not 
been discussed previously in the literature. Before furthering this discussion, we first present in the sequel 
our results. 

In the imaginary-time formalism, the eucledian momentum of the gluon is P^ = {potP) such that 
P^ = (po)^ + P^ with p = p p and pq = 2TrnT where n is an integer. Real-time amplitudes are obtained via 
the analytic continuation po = —iu + 0^ where to is the energy of the gluon. A momentum is said to be soft 
if both Lo and p are of order gT; it is said to be hard if one is or both are of order T. In the Coulomb gauge, 
the complete inverse gluon propagator is: 

D-^{P) = P^S^"' - P^'P•' - ^^"'{P) + —S'''S''^p'p^ , (2) 

where n^''(P) is the gluon self-energy and the last term is due to Coulomb-gauge fixing. We work in the 
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strict 



Coulomb gauge ^c = 0. The gluon self-energy can be decomposed into: 



n^'^(p) = (5n^'^(p) + *ii'"'{p) 



(3) 



where 511 is the hard-thermal- loop and *n is the effective self-energy [4]. P being soft, the hard thermal 
loop is of the same order of magnitude as the inverse free propagator, i.e., (511 ^ (gT)"^, while the effective 
self-energy is of an order of magnitude higher, i.e., *n '^ g{gT)^. Since the momentum running inside 
*n is soft, to calculate it, we have to use effective vertices and propagators instead of their bare (tree) 
counterparts. This ensures the correct expression for the (/(gT)^— correction to the inverse gluon propagator 
and, in particular, that this correction is gauge-independent [4]. 

The hard-thermal- loop SIl is already known from the literature [2,4]. It is real and contributes to the 
determination of the spectrum of the soft gluonic excitations to leading order gT. More explicitly, 611 is 
gauge-invariant and can be expressed in terms of only two independent scalar functions, i.e.,: 



Slf"{P) = SUi{P) ; SU"\P) 

The expressions oi 5IIi{P) and S'nt{P) read [2]: 



-^ m.(P) ; 



p^p^ 



p. 



(4) 



SlliiP) 



Sni{P)^3mlQ,{^) 



SUtiP) = -ml 



fm^\_oJ]Pl\ 



(5) 



\ p J \ p J 

where the Qn are Legendre functions of the second kind. The thermal gluonic mass is given explicitly by 



nig = y/Nc + {l/2)Nf gT/3, where Nf is the number of flavors. 

The effective propagator for soft gluons that intervenes in the calculation of the effective self-energy is 
obtained by inverting (2) while disregarding *n. In the strict Coulomb gauge, its nonzero components are 
*A^°(P) = *Ai{P) and *AJ^(P) = {6'^ -p'p>) *At(P), where *At and *Ai are given by: 

1 . . ,„, 1 



'MP) 



P^-SUt{P) 



'MP) = 



p^-sni{p) 



(6) 



After analytic continuation to real energies, the pole in lu of *At(;) yields the dispersion relation u!t(i){p) of 
the transverse (longitudinal) gluons to order gT. The longitudinal dispersion equation can be rewritten in 
the following explicit manner: 



In ^^^^±^^^(3 
uji- p Sivi 



One finds for soft longitudinal gluons: 



UJI [p) = nig 



1 + lf^V-Af^ 
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i-g/ Z.OU \ni,g/ KJUUU \nLg' 

To order gT, the gluons are not damped. The leading order of the damping rates is obtained after is 
included in the dispersion equations the contribution from the effective self-energy *n, which of course has a 



more complicated structure than that of SIl. One shows that in the strict Coulomb gauge, it also depends on 
two scalar functions such that the mass-shell condition that determines the longitudinal dispersion relation 

is: 

/ _ sni{-ini,p) - *ni{-ini,p) = o , (9) 

where *n; = *n . The longitudinal damping rate is defined by ji{p) = — Im fti{p). Since it is g-times 
smaller than the energy LOi{p), we can write from (9): 

nip) = -Ji — -, — -^ ■ 10 

The denominator in (9) is easy to get from (5). Oncha,s d/dujSIl{—iuj,p)\^^^^^(p\ = —2rag{p/mg)'^-\-. . . . This 
means that in order to get ^i{p = 0), we have to expand the imaginary part of the longitudinal self-energy 
to order p^ . 

In the Coulomb gauge, the only diagrams that contribute to Im*n;(P) above the light cone are the 
three-gluon {ig) and four-gluon [Ag) one loop-diagrams with soft internal momentum [4]. Hence we write: 

Im *ni{P) = - ^ Im Tr,oft [*r°"^"(P, -P, K, -K) * A^J (K) 

+ *T"°^{-Q,P,-K) *A^^'{K) *r^'"'''{-K,P,-Q) *Ay{Q)] , (11) 

where K is the internal loop-momentum, Q = P — K and Tr = PX^fc / (2 )'■* • ^^'^ subscript 'soft' means 
that only soft values of k are allowed in the integral. Eq (11) is what one would normally write for the 3g and 
Ag contributions to the imaginary part of the gluon self-energy, except that everywhere, tree quantities are 
replaced by the corresponding effective ones. The gluonic effective vertices can be written in the following 
form: 

*r(") =r(")-H ,5r(") ; n = 3,4 , (12) 

where the first term is the tree QCD gluon vertex and the second one sums up the contributions from hard 
thermal loops with n external legs. In the case n = 3 it can be written as: 

Jr- {-Q,P,-K)^-Sr'^ i-K^P^-Q) = imgj—^^[—-—), (13) 

where 5* = (i, s) and Vis is the solid angle of the unit three-vector s. Also, PS = ipo + p.s, etc. In the case 
n = 4 we have: 

^ ' ' ' ' aj 4^ psKS \PS-KS PS + KSJ- ^ ' 

To be complete, we give the expression of the 3g tree vertex: 
ri"'^{-Q, P, -K) = -r^'^^(-i^, P, -Q) 



{P + KY S"^ + {Q- KY S^f" -{P + Qf S^"" , 
4 



(15) 



and that of the Ag tree vertex: 



pM^Aff ^p^ _p j^ _^-^ ^ 26'"' S^'^ - (5^^ S"" - S'"' S" 



(16) 



From eq (11) and the discussion before eq (6), we have: 

Im *MP) = -^ lmTj2j^ [*r°"°"(P,-P,X,-X) *MK) 
+ *r"°'^(P, -P, K, -K) {5'^ - fc'F) *At{K) 

+ *r"°"(-Q,p,-i^) *Ai{K) *r"°"(-i^,p,-Q) *A,(g) 

+ *r"°\-Q,P,-K) {6'^ -FF) *At{K) *r^™(-i^,P,-Q) *Ai{Q) 

+ *r-'™(-Q,p,-i^) *AKi^) *r°"*(-i^,p,-Q) (<5*-'' -gV) *At(g) 



+ *r™'(-Q,P,-i^) ((5^^ -FF) *At(ii') *r^'""(-i^,P,-Q) (r5 



^rnn _ A'>T^/^n\ * 



A,(Q)] 



(17) 



There are six contributions: two from the 4g vertices and four from the 3g vertices. We will take henceforth 
nig = 1 and all momenta and energies are in units of it. This simplifies considerably the expressions we 
work with and, if and when needed, the ?ng-dependence can easily be recovered in the final results. Each 
contribution has to be calculated separately. As an illustration, let us discuss how we manipulate the 
contribution from the 3g vertices in which the two effective propagators involved are both longitudinal. The 
corresponding term in (17) is: 



Im *n 



l3gU 



(P) 



g^Nc . „^r- /■ d^fc 



__ *rO°"(-Q,P,-i^) *Ai{K) *r°°"(-X,P,-Q) *A,(Q) 



(18) 



in a clear notation. Using eqs (12), (13) and (15), wc have: 

Im *n,3,„(P) = --g'N^ImTY^J j^ *MK) *A,(Q) 



dfls iko 
~A^PSKS 



dn,^ 



«go 



dris iko 
47r PSKSJ 47r PSQS 



(19) 



We need an expression for the two solid-angle integrals involved in (19). For this purpose, we use the 
simple expansion 1/PS ~ l/ipa (1 — p.s/ipo — p.s^ /p^ + ...), which is valid as long as p < \ipo\, a condition 
satisfied in the region p < nUg before analytic continuation to real energies and after. We then perform the 
angular integrals straightforwardly and get: 



d^s ikn 



ikn 



At: PSKS ipok 

9 



iOk 



px iko 

Vofe — -. — t — r' 

ipo \ k 

- 2; ^ 1 - ^Qok + :^(1 - -^ ) Hr 



iko 
k 



i + %]Qok 



(20) 



where Qofc stands for Qo(^) and x = cos(p,k) . Wc obtain a similar expression for the other solid-angle 
integral: we only need to replace K by Q. 



The next step is to multiply the expressions for the solid- angle integrals as needed in eq (19) while still 
expanding in powers of p. Then we integrate over the solid angle of k and the odd-powered terms in p cancel 
out. For the first piece in (19) we get: 
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(21) 



QoA; 



/vn 



*A;(go,fc) , 
where 9^ stands for d/dk. Similarly, the second piece in (19) reads: 

J (2^ ^'^"^^ ^'^^^ / i^T^^Tf:^ / - 



-9g^NclmTY^ 



47r PSKS J 47r PS'QS' 
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(22) 



3 V ipQ ' ipQ J \ 3 J J k "^2' 
In order to eliminate the appearance of Q^f. in (21) [1], extra work is needed. For this purpose, we use the 
expression for *Ai{K)^^ one obtains using eqs (6) and (5). After some algebra and disregarding terms that 
from the outset do not yield an imaginary part to the effective self-energy after analytic continuation, we 
get: 

-95'7V,Imr^ / j^ *MK) *MQ) 



ka 



(27r) 
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iko 
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(23) 



-(-(- y) (- y (-€))•-'<-'- 5 (' 
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Remember that all momenta and energies are in units of rrig, set to one. Putting the two pieces (22) and 
(23) together, we find: 



Im*n,3<,«(P) 



9.927V, 
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l-2'-^y-^Qokd, (24) 



The next step is to perform the sum over kg. We do this by using the spectral representations [5] of the 
different quantities involved, fn particular, we have for the effective propagators: 



'AmW 



+ 00 



dr e'''" ^ / dcu pt,i (LU,k) (1 + n{uj) ) e~ 



(25) 



where n(cij) = l/(exp(a;/r) — l) is the Bose-Einstein distribution and the spectral densities pt,i{(^, k) have 
the following form: 

Zt.iik) 



Pt,i{uJ,k) = 



S{lo - utAk)) - S{uj + Lut,i{k)) + Pt.A^, k) e(fc2 - Lo^) , 



(26) 



2wt,/(fc) 

where G is the step function. The two residues Zt.i{k) and the two cuts (3t.i{ujTk) are given in eqs (3.23) 
and (3.24) of [1]; see also [5]. The sum over fcg can now be performed straightforwardly. The analytic 
continuation to real energies is obtained by the subsequent replacement ip^ -^ lui (p) + zO^ and we extract 
the imaginary part from our expressions using the known relation l/(x + iO+) = Pr(l/.T) — iirS^x). We find: 



lm*ni3gu{P) = 



3g^N,T 
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dk 



+ OC 
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+ ^ (1 - 2c.i) e,dk + ^ 61 dl 
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In this equation, pu stands for pi{uji,k), i = 1,2, and Oi ~ Q{k^ — cul) . Also, for uj soft, we have used 
n{Lu) ~ T/uj. Finally, note that lui depends on p and in principle eq (27) still needs a further expansion. But 
this will prove unnecessary. 

The other contributions to Im *Ili{P) are calculated in a similar manner. The expressions of the two Ag- 
contributions are quicker to obtain whereas more work is needed for the other 3.g-contributions, especially the 
Sgtt-part. We spare the reader the details and here we give the final expressions. The two 4(7-contributions 
read: 

^n^N.T r°^ r+^ di,u r+°^ di,,^ r 

ujiujik Qi 
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Where e(a;) is the sign function and Si = S{uj'l — fc^). As to the 3g-contributions, we have: 



(28) 



(29) 



Im *ni^gti{P) = Im *ni^gu{P) = 
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and wc have: 
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(31) 



We put now the different contributions together. The terms independent of p caneel out as they should 
and, after reorganizing the terms proportional to p^ in a suitable manner, we obtain, using eq (10) in the 
limit p ^ Q: 



im 



247r 
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where the pure numbers ajp and ajg have the following analytic expressions: 
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Things are now well set up to pursue our initial discussion. When introducing this work, we have 
said that at zero momentum, the longitudinal-gluon damping rate is to be equal to the transverse-gluon 
one, i.e., 7;(0) = 7t(0). If one calculates the transverse-gluon damping rate at zero momentum, one finds 

7t('^) ~ ^24" '^tOi where ato = '^la 1 the expression of which is given in (33) above, see [2]. Therefore, one 

(2) 
automatically expects aj^ to vanish. But this seems not trivial an issue at all. Indeed, take for example the 

term uJi/3 (1 — 2a;i) Oi dkPi2 in the analytic expression of a^Q given in (34). It turns out that the contribution 

(2) 
of this term to a^Q and hence to 7((0) is infrared divergent. More precisely, if one cuts the lower bound of 

the fc-integral to a small momentum ry <C 1, then one straightforwardly works out the contribution of this 



term to ajp and finds: 
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Note that the integral over k in the right-hand side of the above equation is finite. Thus one clearly sees 

that this term is infrared divergent. Other terms have a similar behavior. 

Remember that we decided to carry out this work when we noticed that terms similar to the ones in (34) 

do appear in the analytic expression of a^ in eq (1), terms which are infrared divergent too [1]. At this stage, 

it becomes clear that the appearance of such terms is linked to the expansion in powers of p^. Besides, their 

appearance in the expression of 7;(0) is puzzling to a certain degree because as we said, 7;(0) is expected 

to be divergence-free. It therefore may be that though such terms, when taken individually, are infrared 

divergent, when put together, the infrared pieces cancel out and we are left with a finite contribution, likely 

a vanishing one in the case of a^Q . But it should be clear that checking through this point is a delicate 

matter. In any case, it is beyond the scope of the present work. 

(2) 
If such a cancellation does occur for both a^^' and ati, does it signal to a more basic pattern related 

to some aspects of the infrared behavior of hot QCD or is it a mere 'coincidence' ? On the contrary, if 

the infrared divergences persist, is it possible to devise a sort of finite-temperature infrared- renormalization 

scheme that would absorb them or is it that we will have to 'live with', so to speak? If such a scheme is 

possible, how consistent would it be with other calculations? 

Still in the case where these divergences persist, we can also naturally ask whether they are linked in 

some way to the logarithmic behavior ln(l/g) discussed in [3]; in other words, whether they are merely a 

different manifestation of the same phenomenon. We think it is too early to tackle thoroughly this interesting 

question but nevertheless, we can put forward the following argument. The logarithmic behavior is obtained 

when letting p come from above rrig down into the interior of the soft region, but keeping it larger then what 

is expected to be the magnetic scale, i.e., g^T. Therefore, though arguably we may think of the infrared 

(2) 
divergences in ati as being another expression of this logarithmic behavior, this cannot be argued for oJq 

simply because here p has to be strictly zero and the logarithmic behavior cannot be extrapolated down to 

zero momentum. Besides, Since the infrared divergences we obtain in the case of 7((0) are closer in nature 

to the ones we have for 7t(p), we find it difficult at this stage to answer this question positively, especially 

that the role of the magnetic effects is still somewhat obscure. In any case, one should not forget that the 

quark-gluon plasma is expected to be stable for zero and very small gluonic momenta, and therefore, any 

divergence in the gluonic damping rates, if it persists after scrutiny, has to be adressed from this standpoint. 

It seems to us that all the issues we are raising in this work arc worth investigating, especially that as 
far as we can tell, they haven't been (at least fully) adressed in such a way previously. 

We would like to thank Asmaa Abada and Olivier Pcne from Orsay for interesting and quite useful 
exchanges that actually helped reshape parts of our discussion. 
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